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ABSTRACT 



Aims. We explore how ALMA observations can probe high-redshift galaxies in unprecedented detail. We discuss the main observables 
that are excited by the large-scale starburst, and formulate expectations for the chemistry and the fluxes in the center of active galaxies, 

Oin which chemistry may be driven by the absorption of X-rays. We estimate the expected number of sources at high redshift to infer 
I whether an ALMA deep field may find a reasonable number As a specific example for the complex interpretation of sub-mm line 

V,^ ■ observations, we analyze the recently detected z = 6.42 quasar, for which a number of different line fluxes is already available. We 

note that our diagnostics may also be valuable for future observations in the local universe with space-borne instruments like on 
SPICA or FIRI. 

Methods. To estimate the observables from the starburst, we check which emission from the starburst ring of the nearby Seyfert 2 
^ . galaxy NGC 1068 falls into the ALMA bands if the galaxy were placed at z = 8. We estimate the sizes of the central X-ray dominated 

$H ' region based on a semi-analytic model, and employ a detailed ID approach for the chemistry in X-ray irradiated molecular clouds 

' to evaluate the chemistry and the expected line emission under these conditions. We make use of pre-existing chemistry calculations 

I in X-ray dominated regions to show the dependence of different line fluxes on X-ray luminosity, cloud density and cloud column 

density. We use theoretical models for the high-z black hole population and the local SMBH density to estimate the number of sources 
at higher redshift. 

Results. We show that a number of different fine-structure lines may be used to probe the starburst component of high-redshift quasars 
^ ' in considerable detail, providing specific information on the structure of these galaxies by several independent means. We show that 

oo : the size of the central X-ray dominated region is of the order of a few hundred parsec, and we provide detailed predictions for the 
expected fluxes in CO, [CII] and [01]. While the latter fine-structure lines quickly become optically thick and depend mostly on the 
strength of the X-ray source, the rotational CO lines have a non-trivial dependence on these parameters. We compare our models 
to XDRs observed in NGC 1068 and APM 08279 and find that the observed emission can indeed be explained with these models. 
Depending on the amount of X-ray flux, the CO line intensities may rise continuously up to the (17-16) transition. A measurement 
■ of such high-y lines allows one to distinguish observationally between XDRs and PDRs. For the recently observed z = 6.42 quasar, 

' we show that the collected fluxes cannot be interpreted in terms of a single gas component. We find indications for the presence of a 

dense warm component in active star forming regions and a low-density component in more quiescent areas. Near z = 6, an ALMA 
deep field may find roughly one source per arcmin". At higher redshift, one likely has to rely on other surveys like JWST to find 
appropriate sources. 



Key words. Keywords should be given 

1. Introduction With the upcoming mm/sub-mm telescope ALMAQ, it will 

be possible to probe such structures in even more detail, due to 

The recent detection of kpc-scale star-forming structures its significandy improved sensitivity, angulai- and spectral res- 

at z = 6.42 th rough the detection of [CII] emission olution. From a theoretical point of view, it is therefore inter- 

('Walter et al.','2009a') and emission in various rotational CO lines esting to speculate what ALMA might observe in the center of 

(Riechers et al., 2009) confirms the importance of mm/sub-mm such host galaxies. Within the next ten years, we further expect 

observations to infer gas distribution and dynamics in quasar the advent of SPICA0 and FIRfl which will probe the universe 

host galaxies. Emission in CO and the continuum in high- in the mid- or far-IR regime, respectively. These telescopes will 
redshift quasars hav e also been reported by Omont et al. ( 1996), 

CarilH et al.' ^2002), Walter et al.^2004), Klamer et al. (2005), 

WeiBetal. (2005), Maiolino et al (2007), Walte r et al. (2007) ' http://www.eso.org/sci/facilities/alma/ 

WeiB et al] (|2007|), iRiechers et al.i (l2008bl) and iRiechers et all - http://www.irisas.jaxa.jp/SPlCA/index.html 

t2008a) . ^ http://sci.esa.int/science-e/www/object/index.cfm?fobjectid=40090 
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be able to apply in the local universe what we suggest as high- 
redshift diagnostics for ALMA. 

Before speculating what ALMA may see in the centers of 
high-redshift galaxies, we turn our attention to the properties 
of molecular clouds in the central molecular zone (CMZ) of 
the Milky Way. Studies employing and CO lines indicate 
the presence of high-temperature (T ~ 250 K) and low-density 
(n ~ 100 cm'^ ) gas (lOka et all l2005l) . NH3 observations by 
iNagavama et al.l (l2007l) confirm the presence of warm molec- 
ular clouds, with temperatures mostly between 20 - 80 K. The 
presence of a 120 pc star-forming ring was infered by CO obser- 
vations, indicating typical densities of lO'''''"^ cm"-' and kinetic 
temperatures of 20 - 35 K (iNagai et al.l l2007h . These tempera- 
tures were derived under the conservative assumption of a beam 
filling factor 1, while smaller filling factors would give rise to 
higher temperatures. 

In the centers of active galaxies, the supermassive black 
hole will emit radiation in a broad range of frequencies. 
Particularly interesting is the emission of X-rays, as those pho- 
tons can penetrate molecular clouds even at high column den- 
sities. The resulti ng cloud temp eratures r ange from a hundred 
Kupto 1 000 K (Lepp & Dalgarno, 1996; Mal onev et al.Lll996t 
iMeiierink & Spaans, 2005; Meiierink et al., 200?!) . Such X-ray 
dor ninated regions have b een reported for instance in NGC 1068 
by iQalliano et al.l (l2003l) . Due to the high ambient pressure, 
higher-dens ity clouds m ay form due to the thermal instability 
dWada & Normanil2007h . 

At high redshift, resolution of present-day telescopes is gen- 
erally not sufficient to resolve the central regions of quasar host 
galaxies. Under exceptional circumstances, this is however pos- 
sible. Indeed, highly excited high-7 CO and HCN line emission 
was found in APM 08279+5255 by WeiB et al. (2007.) . As this 
galaxy is gravitationally lensed, it was possible to measure fluxes 
from the central X-ray dominated region that are usually beam- 
diluted. While the CO line fluxes usually rise up to the CO (5-4) 
transition and then decrease, they continue to rise in this system 
up to the CO (10-9) transition, providing clear evidence for the 
presence of warm gas in the central region. Sim ilar results have 
also been obtained for the Cloverleaf quasar ([Bradford et al.L 
I2009h . 

Motivated by these results, we study in more detail the pos- 
sibility to probe high-redshift quasar host galaxies with ALMA, 
and in particular their central regions. In § |2] we review the 
main heating mechanisms that may be present in high-redshift 
quasars, and discuss their influence on the chemistry. In § [3] we 
give a brief summary on the main FDR observables, which are 
already used to probe galaxies at high redshift. As a specific ex- 
ample, we calculate the expected fluxes for the Seyfert 2 galaxy 
NGC 1068 if it were located at z = 8. Our expectations for the 
central X-ray dominated regions (XDRs) are formulated in § |4] 
based on detailed chemical models including ~ 50 species and 
several thousand reactions. Evidence for X-ray dominated re- 
gions at different redshifts is reviewed and discussed in § |5] On 
this basis, we assess the prospects for finding new sources in an 
ALMA deep field in § |6] We conclude in § |7] In summary, this 
paper provides a basic set of predictions concerning observations 
of high-redshift quasars with ALMA. In a companion paper, we 
plan to provide diagnostics based on the observed line fluxes that 
will allow one to infer physical properties such as the star for- 
mation rate or the X-ray luminosity based on the observed line 
emission. 
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Fig. 1. The heating rates per hydrogen atom due to X-ray ab- 
sorption (XDR contribution), absorption of soft UV-photons 
(PDR contribution) and total heating rate as a function of ra- 
dius. The calculation assumes a 10^ Mq black hole, with 3% of 
its Eddington luminosity being emitted in a hard spectral com- 
ponent between 1 and 100 keV with a spectral slope of -1. The 
strength of the soft UV radiation field is taken as Go = 10 in 
Habing units and the effective density «e)f =10^ cm"-'. For such 
a configuration, the XDR contribution clearly dominates within 
the central 100 pc. 
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Fig. 2. The expected size of the X-ray dominated region in pc, 
for a black hole with 10^ Mq with a spectral slope of -1, as a 
function of the soft UV radiation field Go in Habing units and 
the effective number density «efF- The straight lines assume the 
power-law between 1 and 5 keV, while the dotted lines assume 
it between 1 and 100 keV. 



2. Chemistry in high-redshift quasars 

Emission from molecular clouds in active galaxies can be ex- 
cited by a variety of different mechanisms. Mechanical feedback 
may be imp ortant locally, i n part i cular in the presence of sho cks 
or outflows ( Loen en et al 1 12OO8I: IPapadopoulo s et all l2008l) . In 
addition, there is radiation in a broad range of frequencies, both 
from the starburst and the supermassive black hole. UV emis- 
sion generally gives rise to compact HII regions at temperatures 
of ~ 10"* K in which molecules are completely dissociated and 
emission is mostly by Lyman a and various fine-structure lines. 
Such photons are however absorbed by relatively small column 
densities, and stellar HII regions never become larger than a few 
pc, considerably smaller than the scales of interest here. 
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Soft UV-photons have smaller cross sections and may pen- 
etrate larger columns. In the presence of a starburst, they are 
theref ore the dominant driver of t he molecular cloud chem- 
istry (iHollenbach & Tielensl Il999l) . While their heating effi- 
ciency is low (0.1 - 0.3%), they are very efficient in dissoci- 
ating molecules. In such photon-dominated regions (PDRs), one 
generally expects somewhat enhanced cloud temperatures, with 
most of ffie emission in fine-structure lines, as CO is efficiently 
dissociated. A detai l ed rev iew of these processes is given by 
lMeiierink & SpaansI (l2005h . 

X-ray photons have even smaller cross sections than the soft- 
UV photons, and can thus penetrate larger columns. Specifically, 
a I keV photon penetrates a typical column of 2 x 10^^ cm"^, a 
10 keV photon penetrates 4 X 10^^ cm"^ and a 100 keV photon 
9 X 10^° cm"^. For this reason. X-rays can keep molecular clouds 
at high temperatures even at high column densities. They have 
high heating efficiencies of the order of 30%, and are inefficient 
in the dissociation of molecules. A fraction of them may how- 
ever be reprocessed by the gas and converted in soft-UV pho- 
tons, which may lead to some molecular dissociation. Detailed 
reviews are given by Maloney et al. ( 1996); Lepp & Dalgarno 
(Il996h : iMeiierink & SpaansI (|2005j). Therefore, X-ray absorp- 
tion drives a completely different type of chemistry, and may 
potentially result in temperatures up to 1000 K. The fraction of 
molecules, in particular CO, can be very high. 

To estimate the potential extent of such a central X-ray dom- 
inated region, we employ a toy model that takes into account 
the heat input from the starburst and from the X-ray emission 
of the supermassive black hole. We assume here an axisymmet- 
ric situation with a central supermassive black hole and an ex- 
tended molecular disk in the host galaxy. The radiation from 
the SMBH will consist of a soft and a hard component. The 
soft component is easily absorbed at the edge of the molecular 
clouds and can be neglected for column densities of 10^^ cm"^ 
and above. The X-ray photons, on the other hand, penetrate 
deeply into the molecular disk and excite emission there. For 
the X-ray photons, we adopt a power-law spectrum for fre- 
quencies larger than 1 keV and use the cross sections given by 
[Verner & Yakovlev ( 1995). For soft photons from the starburst, 
we adopt a typical frequency of 10 eV and a cross section of 
2.78 X IQ-^^ cm 2 dMeijerink & SpaansLlIOOSh . Typical heating 
efficiencies are 30% in XDRs and 0.3% in PDRs. We assume that 
the radiation field from the starburst is roughly constant within 
the central region. The X-ray radiation field, on the other hand, 
will be geometrically diluted and partially shielded by the gas. 
To calculate the attenuation of X-rays, we introduce an effective 
density for the central region, which is given as 

Heif = alO^ cm-^ + (1 - Q')10^ cm-\ (1) 

where a is the volume-filling factor of dense clouds, 10^ cm"^ 
is a typical cloud density and 10^ cm"^ a typical density of 
the atomic medium. We adopt now a specific reference case for 
which we evaluate the heating rates and the expected size of the 
XDR. The X-ray emission depends essentially on the product of 
black hole mass Mbh, the Eddington ratio Ae, and the fraction 
fx of the total luminosity going into the hard spectrum. We as- 
sume a modest black hole mass Mbh =10^ Mq, an Eddington 
ratio Ae = 30%, which lies in the typical ra nge of measured 
Eddington ratios for high-redshift A. GN (Sha nkar et all l2004t 
iKoUmeier etalL |2006'; 'Sha nkar et all 12009 ). and a fraction of 
fx - 10% of the total luminosity emitted in X-rays. For the 
X-ray spectrum, we adopt a frequency range between 1 and 
100 keV with a spectral slope of -1. 



We further assume that the starburst produces a soft UV- 
radiation field of Go - 10 (Habing units) and an effective density 
of Heff =10^ cm"^^. Such an effective density corresponds to a 
volume-filling factor of order 1 and is therefore an upper limit. 
In real AGN, the effective density may be smaller in the central 
region, implying less attenuation of the X-rays and therefore a 
larger X-ray dominated region. On the other hand, the parameter 
Go may vary as well and depend on the strength of the ongoing 
starburst. For this scenario, the expected heating rates per hydro- 
gen atom are given in Fig.[T] 

To explore the parameter dependence in more detail, we now 
consider the effective density and the strength of the starburst ra- 
diation field as free parameters and check how they influence the 
size of the XDR, keeping black hole mass, Eddington ratio and 
the luminosity fraction in the hard component as specified above. 
We consider two cases, one with the hard component between 1 
and 5 ke V (case A), and one with the hard component between 1 
and 100 keV (case B). The results for the XDR size are given in 
Fig. 121 In case A shielding effects can be clearly recognized and 
the size of the XDR depends more on the effective density than 
on the strength of the starburst. As typical molecular cloud den- 
sities in these environments are ~ 10^ cm"^ and the filling factor 
is of order 1%, we can expect an average density of lO-' cm"^. 
At its largest baseline, ALMA can even resolve spatial scales of 
~ 30 pc at z - 5, and should thus resolve the corresponding 
XDRs. 

Simulations by IWada et aTl ([2009) for the clumpy medium 
at scales o f ~ 30 pc indic a te a v olume-filling factor a ~ 0.03. 
Models by iGaUiano etal] (l200l for NGC 1068, on scales of 
a few hundred parsec, indicate a volume-filling factor of 0.01, 
which still leads to a surface-filling factor of order 1. For the 
central 200 p c of our galaxy, values of a ~ 0.1 - .01 have been 
suggested (.McCall et all 119991: lOka et all l2005h . For a large 
volume-filUng factor, we expect a somewhat smaller XDR due to 
the attenuation of X-rays. At the same time, however, this region 
will consist of a large number of clouds that are highly excited. 
For smaller clumping factors, the number densitiy of clouds will 
be reduced, but the size of the XDR increased. 

Of course, the approach used here is only approximate, as 
the attenuation may depend on the direction and be stronger in 
some directions and weaker in others. However, these order-of- 
magnitude estimates should still be applicable for a broad range 
of conditions and also hold in cases of spherical rather than flat- 
tened structures. An implicit assumption of the model is that we 
average over sufficiently large scales where the effective den- 
sity provides a good approximation with respect to X-ray atten- 
uation. For more clumpy structures, the XDR would be more 
inhomogeneous and reach out further along the low-density re- 
gions. In case the large-scale dumpiness is considerable, so that 
molecular clouds do no longer fill the projected surface area in 
the beam, our model predictions need to be corrected with the 
corresponding area-filling factor. 

3. Observables in the PDR 

In photon-dominated regions, soft UV-photons from the star- 
burst provide some heat input for molecular clouds, in partic- 
ular at low column densities, and are efficient in dissociating 
molecules like H2 or CO. The main coolants in this regime 
are therefore fine-structure lines of [CII] and [OI]. Depending 
on the strength of the radiation field, some flux may also be 
emitted in molecules like CO or HCN, in particular in the low- 
lying rotational transitions, and there are further fine-structure 
lines that may contribute as well. As mentioned in the intro- 
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duction, there are plenty of observations that studied emission 
from starburst galaxies with present-day sub-mm telescopes 
(e.g. Omont et al., 1996; Carilli et al., 2002; Walter et al., 2004; 
Klameretal., 2005; WeiB et al., 2005; Maiolino et al., 2007; 
Walter et al.^ ,2007i: iRiechers et ali i2008b.a; ,Greve et al.. ,2009c 
Riechers et al.L 120091; IWalter et al.Ll2009al) 



Observable 



A [/jm] (fi [mjy] min. Redshift 



[O I] 'Pi P, 


63.2 


-0.18 


z>5.7 


[O III] Po 


51.8 


-0.18 


z>3.8 


[N II] ^2 -^^ Pi 


121.9 


-0.07 


z > 2.5 


[O I] ^3 Pi 


145.5 


-0.03 


z> 1.9 


[C II] 'P3,2 PU2 


157.7 


-0.7 


z> 1.7 


[S III] V, ^3 


33.5 


- 1.5 


z > 11.1 


[Si II] '-P3,2 Pi/2 


34.8 


-0.06 


z> 11.1 



quency bands, and what would be the exp ected amount of flux. 
The fluxes given bv'So inoglio et alJ (l2005l) have been measured 
with a frequency resolution of 1500 km/s. Correspondingly high 
velocities can indeed be reached in the presence of fast jets or 
outflows. However, typical line profiles show that most of the 
flux is in a range of +150 km/s, which we adopt here as a fidu- 
cial value. The so obtained observable line transitions are sum- 
marized in Table [1] while the expected sensitivity and angular 
resolution at the largest baselines is given in Table |2] For all the 
transitions, a 3cr detection seems possible for an integration time 
of a few hours. The table does not include CO transitions, as the 
PDR would only excite the low-lying transitions which would 
not fall in ALMA's frequency range for z > 8. We note that the 
expected ratio between [Nil] and [CII] is comparable t o the ob- 
servational upper limit derived bv lWalter et al.l (l2009bl) . 



Table 1. The main observables for ALMA in PDRs at high red- 
shift. This specific example assumes a galaxy with a starburst as 
in NGC 1068 placed at z = 8. We also give the minimal redshift 
from which the lines would be redshifted into the ALMA bands. 



With the upcoming sensitivity of ALMA, we expect that 
such PDRs can be probed in more detail as well. Therefore, 
bright lines like the [CII] 158 fim line can be detected at higher 
significance, allowing higher spectral resolution and probing the 
velocity structure of the gas in more detail. Also, weaker lines 
may be detected as well, probing gas at different densities and 
providing additional information on the chemical conditions. 



4. Expectations for the XDR 

As ALMA may for the first time detect and resolve emission 
for the central X-ray dominated regions, we want to assess here 
in more detail the expected chemical conditions in the central 
region and the corresponding fluxes in different lines. We start 
by discussing the implications of X-rays for the conditions in 
molecular clouds. We then show how ALMA observations can 
distinguish between X-ray chemistry and an intense burst of star 
formation on the same spatial scales. Afterwards, we provide a 
set of systematic model predictions, first assuming an XDR of 
constant size, but also considering the potential increase of the 
XDR in case of a higher X-ray luminosity. 



band 


freq. [GHz] 


eres ["] 


Sc [miy] 


Si [mJy] 


^'beara ["] 


3 


84-116 


0.034 


0.019 


0.163 


56 


4 


125 - 169 


0.023 


0.023 


0.174 


48 


5 


163-211 


0.018 


0.298 


2.63 


35 


6 


211-275 


0.014 


0.039 


0.225 


27 


7 


275 - 373 


0.011 


0.077 


0.372 


18 


8 


385 - 500 


0.008 


0.143 


0.620 


12 


9 


602 - 720 


0.005 


0.232 


0.813 


9 



Table 2. Frequency range, angular resolution ^res at the largest 
baseline, line sensitivity S / for a linewidth of 300 km/s and con- 
tinuum sensitivity Sc for 3cr detection in one hour of integration 
time and primary beam size 0beam- 3 more bands might be added 
in the future, band 1 around 40 GHz, band 2 around 80 GHz and 
band 10 around 920 GHz, which will have similar properties as 
the neighbouring bands. 



To obtain a rough estimate on the expected PDR fluxes in 
different lines, we have evaluated the PDR fluxes that we would 
expect for a system like the Seyfert 2 galaxy NGC 1068, if it 
were placed at high redshift. We adopt z = 8. This system con- 
sists of a central X-ray dominated region dGalliano et al.L l2003h 
and a circumnuclear starburst ring of ~ 3 kpc in size, with a stel- 
lar mass of ~ 10^ M© and an age of 5 Myr (*Spinogl io et al.L 
[2005). On scales of a few hundred parsecs, on e finds a star 
formation rate of a few times 10 Mq yr"' kpc"^ dDavies et al.L 
l2007i) . This is close to the star formation rate in Eddington- 
limited starbursts as suggested by [Thompson et al.l (l2005h . To 
understand which of the fluxes emitted in this region would be 
detectable with ALMA if this system we re located at z 8, we 
went through the spectroscopic sample of lSpinogho et al. I (f200l 
and checked which lines would be redshifted into the ALMA fre- 



4. 1 . Implications of X-rays for molecular clouds 

We follow the chemistry in a one-dimensional molecular 
cloud complex irradiated by X-rays with the XDR code of 
iMeiierink & SpaansI (120051) . The model includes more than 
50 chemical species and several thousand reactions. For CO, 
the detailed level populations are solved consistently with 
the ID radiation transport equation (iPoehnan & SpaansL I2005L 
2006). As the low-metalli city case was explored in detail by 
Spa ans & MeiierinkI (l2008h . we focus here in particular on sit- 
uations with about solar metallicity. 

The first model we discuss corresponds to the XDR in the 
Seyfert 2 galaxy NGC 1068 (see § 15. lb . We plot chemical abun- 
dances and CO emission for cloud column densities between 
10^" - 10^"* cm"^ in Fig. [3] Larger column densities correspond 
to extreme ULIRGs like Arp 220 that is even optically thick 
around 350 GHz (P. Papadopoulos, private communication). The 
fiducial gas density of 10^ cm"-* has little impact on our re- 
sults, unless it drops to below lO"*'^ cm""*. Above this limit, 
the XDR properties are determined by the ratio of X-ray flux 
to gas density. For lower densities, emission in the high-7 CO 
lines would not be excited due to the critical densities. However, 
high-density gas app ears to exist in the center of NGC 1068 
dGalliano et al.Ll2003h . 

The strong X-ray flux of ~ 170 erg s ' cm ^ in NGC 1068 
suffices to make the gas essentially atomic and leads to high tem- 
peratures of ~ 3000 K, as well as relatively low CO abundances 
of the order 10"^ due to photodissociation by soft UV-photons 
produced after the absorption of X-rays. However, the CO abun- 
dance is still higher than in typical PDRs, and the CO intensity 
is high, due to the strong thermal excitation in the hot gas. For a 
column of 10^^ cm ~^, our results appear o f the same magnitude 
as in the model of iGalhano et al.l d2003h . For larger columns, 
the temperature gradually decreases, the gas becomes molecular 
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Fig. 3. A model for the X-ray chemistry in NGC 1068. The 
adopted flux impinging on the cloud is 170 erg s"' cm"^. The 
adopted density is 10^ cm"^. Top: The abundances of different 
species as a function of column density. Middle: The low-J CO 
lines as a function of column density. Bottom: The high-7 CO 
lines as a function of column density. 



and CO gets more abundant, and we find intensities of the order 
10"^ erg s ' cm"^ sr ' in the high-/ CO lines. 

To explore the dependence of the chemistry on the X-ray 
flux, we consider two additional cases. An extreme case with 
~ 1 erg s ' cm"^ is shown in Fig.H] In this model, we find lower 
temperatures of ~ 70 K, a large fraction of molecular gas and CO 
abundances of the order of 10 While the lower temperature 
tends to decrease the CO line intensities, they are still enhanced 
due to the larger CO abundance. Above a column of 10^^ cm"^, 
the intensities increase rather slowly as the lines become opti- 
cally thick. 

As an intermediate scenario, we consider a source with an 
X-ray flux of ~ 10 erg s"' cm"^ (see Fig. |5]l. In this model. 




20 21 22 23 24 

Logio (Column density [cm~^]) 



Fig. 4. The X-ray chemistry in a system with X-ray flux of 
1 erg s ' cm"^ impinging on the cloud. The adopted density is 
10^ cm"^. Top: The abundances of different species as a function 
of column density. Middle: The low-7 CO lines as a function of 
column density. Bottom: The high-/ CO lines as a function of 
column density. 

the temperature is increased to ~ 100 K. The CO abundance is 
initially of the order 3 x 10 and increases to ~ 10 for larger 
columns. For columns less than 10^^ cm"^, the intensities are 
thus reduced by about an order of magnitude compared to the 
previous case, while they are increased by an order of magnitude 
for larger columns. 

4.2. Separating the XDR from a nuclear starburst 

In the center of an active galaxy, not only the X-ray emission 
is enhanced, but one may expect the presence of a strong nu- 
clear starburst. For instance, Arp 220 harbors such a starburst 
on scales of ~ 300 pc. We therefore compare the expected CO 
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Log,() (Column density [cm ^]) 

Fig. 5. The X-ray chemistry in a system with X-ray flux of 
10 erg s~' cm~^ impinging on the cloud. The adopted density is 
10^ cm"^. Top: The abundances of diff'erent species as a function 
of column density. Middle: The low-7 CO lines as a function of 
column density. Bottom: The high-7 CO lines as a function of 
column density. 

line SED of a strong starburst with Go - 10^ with the CO line 
SEDs in X-ray dominated regions, based on the models provided 
bv .Meiierink et alj (l2007h . We normalize them such that the CO 
(10-9) transition has the same intensity in all models. In this 
case, the SEDs can hardly be distinguished at the low-7 tran- 
sitions that are typically observed at low redshift (see Fig.|6|l. At 
higher-7 transitions, the PDR SED drops considerably and flat- 
tens on a low level due to the small amount of hot gas in the 
outer layer of the molecular cloud. We expect that a value of 
G() - 10^ is a robust upper limit for the soft-UV flux that can be 
obtained in a galaxy. In fact, larger values have never been in- 
dicated in previous observations, and indeed such a value would 
require extreme conditions as in the Orion Bar throughout all 
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Fig. 6. A comparison of the CO line SED in case of an intense 
starburst with Go = 10^ with the corresponding SED for molec- 
ular clouds under X-ray irradiation, for different X-ray fluxes in 
erg s"' cm"^. The spectra are normalized such that they have 
the same intensity in the 10th transition. If the impinging X-ray 
flux is at least 2.8 erg s"' cm"^, observations of the 15th CO 
rotational transition can clearly discriminate between PDR and 
XDR chemistry. 



of the galaxy. In XDRs, a much larger fraction of the gas is at 
high temperatures, and thus the SED is not expected to drop as 
rapidly. 

If the total amount of energy injected by X-rays and by soft- 
UV photons is comparable (within a factor of 10), then the pres- 
ence of X-rays can be clearly inferred using the CO (16-15) 
transition, if the local X-ray flux is at least 2.8 erg s"' cm"^. 
However, as discussed in § |2] we expect X-ray flux to dominate 
over the soft-UV in the center of the galaxy. Observations at even 
higher-y transitions may be useful to determine the local amount 
of X-ray flux from the CO line SED. A potential uncertainty is 
the presence of cold dust, which may to some degree absorb the 
CO line emission and thus change the appearance of the SED. 
Due to the characteristic scaling of dust absorption with wave- 
length, we expect that such a behavior could be recognized and 
potentially corrected. For this purpose, it is of course desirable 
to measure as many high-7 CO lines as possible. 

4.3. Model predictions for XDRs of constant size 

Although we have shown in §|2]that the central XDRs can likely 
be resolved with ALMA, it is currently not clear how the ex- 
pected XDR size varies with X-ray luminosity. If the strengths 
of the soft-UV field is independent of this, one should on average 
expect a larger XDR for higher X-ray luminosities. However, it 
is also conceivable that the X-ray luminosity is indicative of the 
system as a whole, and that a higher X-ray luminosity may be 
accompanied by a stronger soft-UV field. In such a case, one 
might expect a smaller increase in the XDR size or even a con- 
stant size. For this reason, we will consider two extreme cases, 
assuming that more realistic scenarios should lie in between the 
two. In this subsection, we will assume that the size of the XDR 
is always constant, of ~ 200 pc. Of course, the numbers given 
here can be easily rescaled for other XDR sizes, or for area- 
filling factors smaller than 1. In the following subsection, we 
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Fig. 7. The expected flux in mJy for high-/ CO lines, for a cen- 
tral XDR of 200 pc, and molecular clouds of 10^ cm"^, as a func- 
tion of X-ray luminosity and cloud column density. We focus on 
lines that fall in ALMA band 6, which offers a good compromise 
between angular resolution and sensitivity. For a source at z = 5 
(soUd line), this corresponds to the (10-9) CO transition, for a 
source at z = 8, it corresponds to the (14-13) CO transition. 



Fig. 9. The expected flux in mJy for the [CII] 158 /im line, for 
a central XDR of 200 pc, and molecular clouds of IC* cm"^, as 
a function of X-ray luminosity and cloud column density. The 
solid line corresponds to a source at z = 5, and the dashed line 
to a source at z = 8. At z = 5, the line is redshifted into ALMA 
band 7, with a sensitivity of 0.08 mJy (1 day, 3cr, 300 km/s). At 
z = 8, it falls into ALMA band 6 with a sensitivity of 0.04 mJy. 




Fig. 8. The expected flux in mJy for the [CII] 158 fim line, for 
a central XDR of 200 pc, and molecular clouds of 10^ cm"^, as 
a function of X-ray luminosity and cloud column density. The 
solid line corresponds to a source at z = 5, and the dashed line 
to a source at z = 8. At z = 5, the line is redshifted into ALMA 
band 7, with a sensitivity of 0.08 mJy (1 day, 3cr, 300 km/s). At 
z = 8, it falls into ALMA band 6 with a sensitivity of 0.04 mJy. 



Fig. 10. The expected flux in mJy for the [OI] 63 //m line, for 
a central XDR of 200 pc, and molecular clouds of 10^ cm"^, as 
a function of X-ray luminosity and cloud column density. The 
solid line corresponds to a source at z = 5, and the dashed line 
to a source at z = 8. At z = 5, the line is redshifted into ALMA 
band 10, with a sensitivity of 0.2 mJy (1 day, 3cr, 300 km/s). At 
z = 8, it falls into ALMA band 8 with a sensitivity of 0.13 mJy. 



will then discuss the implications of varying the X-ray luminos- 
ity for constant Gq. 

For a first estimate, let us consider at source at z = 5 with 
an XDR of at least 100 pc, corresponding to an angular scale 
of 0.016", and a typical intensity in the high-/ CO lines of 
lO""* erg s"' cm"^ sr"'. As can be seen from the calculations 
above, such an intensity can be reached in a broad range of sys- 
tems for column densities of at least 10^^ cm"^, and in fact also 
for column densities of at least 10^^ cm"^ in the presence of suffi- 
cient X-ray flux. With a fiducial velocity dispersion of 300 km/s, 
this corresponds to a flux of 0.03 mJy, which is detectable in a 
bit more than a day in ALMA band 6. 



In a similar way, it is possible to calculate the expected fluxes 
also in various fine-structure lines. Based on the detailed param- 
eter study provided by Meiierink et alj (|2007|P I which shows the 
expected fluxes in various lines as a function of X-ray luminosity 
density and column density, we therefore provide detailed pre- 
dictions for the fluxes from the central X-ray dominated region, 
assuming a characteristic size of 200 pc, consistent with the re- 
sults obtained in § |2] 

For the high-/ CO lines, we focus on those which are red- 
shifted into ALMA band 6, which offers a good compromise 
between angular resolution (0.014" at the largest baseline) and 

http://www.strw.leidenuniv.nl/ meijerin/grid/ 
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Fig. 11. The expected flux in mJy for the [OI] 63 /jm line, for 
a central XDR of 200 pc, and molecular clouds of 10"* cm"^, as 
a function of X-ray luminosity and cloud column density. The 
solid line corresponds to a source at z = 5, and the dashed line 
to a source at z = 8. At z = 5, the line is redshifted into ALMA 
band 10, with a sensitivity of 0.2 mJy (1 day, 3cr, 300 km/s). At 
z = 8, it falls into ALMA band 8 with a sensitivity of 0.13 mJy. 
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Fig. 12. The expected flux in mJy for the [OI] 146 /im line, for 
a central XDR of 200 pc, and molecular clouds of 10^' cm"^, as 
a function of X-ray luminosity and cloud column density. The 
solid line corresponds to a source at z = 5, and the dashed line 
to a source at z = 8. At z = 5, the line is redshifted into ALMA 
band 7, with a sensitivity of 0.37 mJy (1 day, 3o-, 300 km/s). At 
z = 8, it falls into ALMA band 6 with a sensitivity of 0.225 mJy. 



sensitivity (~ 0.04 mJy for an integration time of one day, 3cr 
detection and a line width of 300 km/s). At a redshift z = 5, 
this corresponds to the (10-9) CO transition, at z = 8 to the (14- 
13) CO transition. We assume that the projected surface area is 
homogeneously filled with molecular clouds with central den- 
sities of 10^ cm"^. Fig. |7] shows how the expected flux varies 
as a function of the cloud column density and the X-ray lumi- 
nositjO. Similar results would be obtained for cloud densities of 



^ For the conversion between X-ray luminosity and flux, we assumed 
optically thin conditions. One may need to correct for further attenua- 
tion in case of large filling-factors. 



10"^'^ cm""*, while it is more difficult to excite CO emission at 
lower densities. 

The figure illustrates that higher fluxes can be obtained for 
larger column densities, while intermediate X-ray fluxes are 
ideal for stimulating emission in the CO lines considered here. 
This is because at very high fluxes, a significant amount of X- 
rays would be converted into soft UV-photons and dissociate the 
molecules. 

In Figs. [8] and |9] we show the corresponding results for the 
[CII] 158 line, both for a density of 10^ cm""* and a density 
of 10^ cm"^. For densities of 10^ cm"-^, low column densities 
are sufficient to yield a detectable amount of flux, and the flux 
monotonically increases with the X-ray luminosity. As in the 
PDR case, this line is therefore valuable to explore the centers 
of high-redshift quasars, and provides complementary informa- 
tion to the CO lines. For densities of lO'* cm"^, we find a stronger 
dependence on column density, and indeed columns of at least 
10^^ cm"^ are needed to yield a detectable amount of flux. 

The results for the [OI] 63 /im line are given in Figs. [TOl and 
[TTl again for densities of 10^ cm"-' and 10^ cm"-', respectively. 
This line quickly becomes optically thick. Therefore, in particu- 
lar for cloud densities of 10^ cm"^, it is insensitive to the column 
density, but provides a good measure for the X-ray flux. It is very 
bright. Even for cloud densities of 10"^ cm"^, it depends more on 
X-ray luminosity than column density, and is still detectable in 
the case of high column densities and strong X-ray fluxes. It is 
thus well-suited to study gas dynamics in the central XDR by 
resolving the line profile. 

We also provide results for the [OI] 146 yum line in Fig. [12] 
for a cloud density of lO*" cm"^. For lower densities, this line 
is hard to excite. It can also be very bright and shows a strong 
dependence on the X-ray flux. The ratio between the [OI] 63 jjm 
line and the 145 //m line is generally about 0.1. 

Emission from neutral carbon seems more difficult to detect. 
The intensity of the [CI] 369 /im line is typically at least an order 
of magnitude smaller than the intensity in the [CII] 158 //m line, 
and only significant in the case of strong X-ray fluxes and high 
column densities. As discussed in § 15.31 the relatively low ra- 
tio between these lines in the recently detected z = 6.42 quasar 
puts a significant constraint on theoretical models. Other car- 
bon lines like [CI] 609 /im or [CI] 230 /zm are even weaker and 
should never be visible. Additional lines like [Sill ] 35 or 
[Fell] 26 //m can be bright as well (Meiie rink et all 12007 ). but 
typically have too short wavelengths for ALMA, except at red- 
shifts z > 9. 

The reader may notice that for the model predictions given 
in this subsection, we restricted ourselves to a range of X-ray 
luminosities of about two orders of magnitude, c orresponding to 
the ran ge of X-ray fluxes avaiable in the data bv lMeiierink et al.l 
(I2007h . This range of data has been chosen such that it covers 
the observationally interesting cases. For lower fluxes, or in our 
case lower X-ray luminosities, we would not expect significant 
emission driven by X-rays, as can be seen in the corresponding 
figures. Similarly, it is straightforward to extrapolate the behav- 
ior towards higher X-ray fluxes: As shown in Fig. |2l the CO 
intensities decrease considerably at higher luminosities. This is 
because part of the X-rays are reprocessed to soft-UV photons 
that dissociate all the CO. At even higher intensities, the gas tem- 
perature would increase up to lO'* K and become fully ionized. 
Molecules would no longer survive under such circumstances. 
The emission of [CII] may increase a bit further, until a large 
fraction of carbon is doubly-ionized at temperatures near 
Similar considerations hold for the oxygen line. However, it is 
not clear whether such scenarios are actually physical, or if the 
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Fig. 13. The expected flux in mJy for the CO (14-13) transition 
(solid line), [CII] 158 fim (dotted line), [OI] 63 /im (dashed line) 
and [OI] 146 /urn (dot-dashed line) emission for a quasar at z = 
5, as a function of X-ray luminosity and average density. We 
assume a typical cloud column density of 10^^ cm"-, and a soft- 
UV field Go = 100. 



Fig. 14. The expected flux in mJy for the CO (14-13) transition 
(solid line), [CII] 158 fim (dotted line), [OI] 63 /urn (dashed line) 
and [OI] 146 /^m (dot-dashed line) emission for a quasar at z = 8, 
as a function of X-ray luminosity and cloud column density. We 
assume a typical cloud density of lO'* cm"^, and a soft-UV field 
Go = 100. 



molecular cloud would rather be photo-evaporated at that point. 
As already mentioned above, the size of the XDR may be more 
extended in the case of such high luminosities. Then, moderate 
X-ray fluxes will be present on larger scales. 

In summary, we can therefore conclude that the main observ- 
ables for the central XDR are the high-7 CO lines as well as the 
fine-structure lines of [CII] and [OI]. The [OI] 63 and 146 jjm 
lines show a strong dependence on X-ray luminosity and may 
provide a good handle on this quantity. In combination with an 
observation of the [CII] or high-7 CO lin es, the X-ray luminosit y 
can be estimated as well. As shown by iMeijerink et al.J (|2007|), 
such line ratios also provide valuable information on gas den- 
sity and temperature. We also note that it is possible to discrim- 
inate such XDRs from reg ions with strong mechanical heating 
dPapadopoulos et al 1 12008). This can be done for instance on 
the basis of the observed X-ray luminosity, or by looking at the 
dust SED. While in XDRs, both dust and gas will be at high 
temperatures in dense clouds, there should be a clear discrep- 
ancy between gas and dust temperature if heating is due to local 
shocks. We therefore expect that the physical conditions in the 
central XDRs can be probed in detail with ALMA. 

4.4. Model predictions for variable XDR-sizes 

As mentioned above, the size of the XDR may increase consider- 
ably for increasing X-ray luminosity. We first explore the depen- 
dence on the average density and the X-ray luminosity, assuming 
a typical cloud column density of 10^^ cm"^, and a soft-UV field 
Go - 100. The X-ray spectrum is assumed to range from 1 to 
100 keV. As shown in §|2] the size of the XDR then varies as a 
power-law with density. The same is true for the expected fluxes 
of CO, [CII] and [OI], as shown in Fig. [13] We estimate those by 
adopting a typical value in mJy on scales corresponding to the 
size of the XDR. These should indeed consitute the main contri- 
bution. For [CII] and [OI], the total flux may be a bit larger, as 
the expected amount of emission increases with the X-ray flux 
(see previous subsection). CO, on the other hand, might be ef- 
ficiently dissociated in the very inner core. However, we do not 
expect this to affect our predictions significantly. 



We also explore the role of the molecular cloud column den- 
sity for the expected fluxes. For this purpose, we adopt an av- 
erage density of lO'* cm"^ and a soft-UV field Go = 100. In 
Fig. [14] we show the expected fluxes as a function of X-ray lu- 
minosity and column density. As a generic feature, we find that 
the expected fluxes vary for columns smaller than 10^^ cm"^, but 
level-ofF for higher values. The fluxes in high-column density 
systems should thus essentially depend on the X-ray luminosity 
only. 

5. Evidence for XDRs and the interpretation of 
sub-mm line observations 

Although the presence of a central X-ray dominated region 
seems unavoidable from a theoretical point of view, we want 
to review current evidence for the presence of central XDRs 
in active galaxies. Such evidence is present in local galaxies 
like NGC 1068 that can be studied in great detail, as well as 
in high-redshift quasars like APM 08279 in which the corre- 
sponding fluxes are magnified by a gravitational lens. Similar 
indications are present also in the Cloverleaf quasar, where CO 
fluxes up to the (9-8) transition have been m easured, and no 
turnover in the line SED has been found yet ([Bradford et all 
|2009). These sources will allow a first test for the expectations 
we have formulated in the previous section once the ALMA tele- 
scope becomes available. We also discuss the z - 6.42 quasar 
SDSS Jl 14816.64+525150.3, as it provides a good and interest- 
ing example concerning the complex interpretation of sub-mm 
line observations. We will further discuss how its properties can 
be understood better if additional data are provided, with partic- 
ular focus on the role of ALMA. 

To distinguish between different excitation mechanisms of 
molecular clouds, it is very important to have observational 
diagnostics for the various excitation mechanisms. First ef- 
forts for modelling the che mistry in XDRs w e re per formed 
by Malo nevetaP j 19961) and Lepp & Dalgarno (1996), while 
PDR c hemistry was originally studied by Hollenbach & Tieleni 
(Il999h . Recent effort s to discriminate s uch models have 
been performed by iPerez-Beaupuits et al.l (l2007l l2009l) . and 
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Fig. 15. A sketch for a situation with an inhomogen e ous X DR, 
motivated by the observations of iGalliano et alj (l2003l) in 
NGC 1068. While X-rays are shielded by a central absorber 
along the line of sight, they may stimulate emission in molecular 
clouds in the perpendicular direction with the typical character- 
istics of XDRs. 



new diagnostic diagrams to discriminate AGN- and starburst- 
do minated ga laxies have been provided by ISpoon et"an (l2007h 
and lHao et al.i (,2009.) . 



5.1. NGC 1068 



As shown bv lGalhano et al.l(l2003h . NGC 1068 contains a central 
XDR. Their analysis is based on the observ ed intensitiy of the 
H2 2.12 f jm line ( Galliano & Alloin, l2002h and the rotational 
CO lines jSchinnerer et al.l, i2000). and showed that the observed 
emission can be cons istently explained with the XDR-model of 
iMalonev et al.l (Il996h under the following assumptions: 

- The central engine is a power-law X-ray source with spec- 
tral slope a - -0.7 and luminosity of 10"''* erg s"' in 
the 1-100 keV range, consistent with the X-ray lumi- 
nosity determined fro m VLBl water maser observations 
(iGreenhill et al.L [T996'). 

- The emission originates from molecular clouds with a den- 
sity of 10^ cm"-', a column of 10^^ cm"^ at a distance of 70pc 
and solar metallicity. 

They show that the central XDR is indeed inhomogeneous 
due to a central X-ray absorber that shields the X-rays along the 
line of sight, while they can stimulate molecular emission in the 
perpendicular direction. A sketch of such a situation is given in 
Fig. [15] In the presence of a torus, indeed it seems likely that X- 
ray emission is preferred in the direction orthogonal to the torus, 
and the XDR fluxes may therefore be preferentially detectable in 
situations where the molecular disk is observed face-on. 



5.2. APM 08279 

In the z - 3.9 galaxy APM 08279, the usual geometric effects 
concerning emission from the central regions are compensated 
by a gravitational lens. It therefore provides an ideal test case 
to study the emission we might see when future telescopes like 
ALMA can observe the central regions with higher sensitivity 
and higher spatial resolution. As expected for molecular clouds 
excited by X-rays, they find that the CO line intensity increases 
up to the (10-9) transition. Significant flux is also detected in the 
HCN(5-4) line. 

Based on brightness temperature arguments, the results from 
high-resolution mapping and lens models from the literature. 



IWeiB et al.l (l2007h show that the molecular lines arise from a 
compact (100-300 pc), highly gravitationally magnified (m=60- 
110) region surrounding the central AGN. It is interesting to note 
that this amount of magnification is comparable to the increase 
in sensitivity due to ALMA. They can distinguish two compo- 
nents of the gas: A cool component with a density of 10^ cm"^ 
and a temperature of 65 K, on scales larger than 100 pc, and 
a warm component of gas with lO"* cm"-' and temperatures of 
220 K on scales smaller than 100 pc. So far, their results do not 
provide indications for inhomogeneities. 

We have c ompared the turn over for the rotational CO 
lines found by IWeifi et all (12007 ^ with the grid of PDR and 
XD R model predictio ns that were made publicly available by 
Me iierink et al.1 (l2007h . In the PDR case, we find that cloud den- 
sities of 10^ cm"^ and a radiation field of Go ~ lO"*^^ with a 
cloud column of ~ 6 x 10^^ cm"^ are required. For the XDR 
case, the turn-over can be explained with a cloud density of 
lO"*'^^ cm""*, an X-ray flux of 2.8 erg s"' cm"^ and a cloud col- 
umn density of ~ 2.6 x 10^^ cm. In both cases, we need to require 
that additional cold gas is present that emits in particular in the 
low-lying J levels to explain the observed ratio between high-7 
and low-7 rotational CO lines. 

The X-ray spectrum of AP M 08279 has been observed with 
Chandra dChartas et al.L l2002h . indicating a luminosity of 4 x 
lO'^^m"' erg s"'. Adopting a lens magnification of m ~ 100, one 
expects a flux of ~ 354 erg s"' cm"^ at a distance of 100 pc for 
an optically thin source. However, for a density of lO"*'^^ cm"^, 
we expect significant attenuation effects th at may considerabl y 
decrease the ffux in the cloud. As shown bv lWada et al.l ( l2009l) . 
shielding may locally vary by two orders of magnitude due to 
column density ffuctuations in the torus and give rise to a peak 
optical depth of t ~ 5 for frequencies of 3 keV. 

With ALMA, it will be possible to probe the distribution of 
these gas components and their kinematics in even more detail. 
Due to the higher sensitivity, the error bars in the flux mea- 
surement will decrease further and one may probe whether the 
turnover occurs at the (10-9) transition, as indicated now, or in- 
deed even at higher-7 levels. In this case, one could clearly dis- 
criminate between PDR and XDR models. 

5.3. SDSSJ1 1481 6.64+5251 50.3 

The discovery of the extended^ source 

SDSS Jl 14816.64+525150.3 at z = 6.42 by IWalter et al.1 
(I2009ah has stimulated great observational intere st that led to a 
broad collection of data on this source. Recentlv. iRiechers et al.l 
(2009) provided a table with fluxes in the CO (3-2), (6-5) and 
(7-6) transitions, fluxes in the CI line and the [CII] line, as well 
as upper limits on five additional lines. Although the publication 
of such upper limits is very valuable and can provide important 
constraints on theoretical models, for simplicity we focus here 
on those lines that were actually detected. As we shall see, the 
interpretation of those is akeady complex, and a more thorough 
analysis would be beyond the scope of this work. The models 
and data emp oyed here are b ased on the publicly available grids 
provided by Meiierink et al.l (l2007i) . Due to the large scales of 
this source, we assume that most of the flux is driven by PDR 
chemistry, especially because the relevant scales for XDRs are 
currently unresolved. 

To understand the complexity of these data, it is illustrative 
to check whether a model with one given density and a fixed pa- 
rameter Go is able to explain them. For CO, the ratio between 
the (6-5) and the (3-2) transition is about 3.35, while the ratio 
between (7-6) and (3-2) is about 3.15. Although it appears that 
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the turnover in the SED may have been reached, this is not fully 
clear from the data, as the uncertainty in the fluxes is ~ 10%. 
The increased intensity in the (6-5) transition requires a signifi- 
cant amount of soft-UV flux, while the almost flat behavior indi- 
cates that the gas should be at intermediate densities of about 
~ lO'* cm"^. At higher densities, the seventh rotational level 
would be more easily excited, which would give rise to a sig- 
nificant discrepancy between the (6-5) and the (7-6) transition, 
even for low values of Go- For lower densities, on the other hand, 
it becomes diflicult to excite these rotational levels. A reasonable 
fit to the CO SED can thus be obtained with n - lO'* cm""* and 
Go = 10'". 

A problem arises, though, if the observed fluxes in [CII] and 
CI should arise from the same gas component. In this case, our 
model predicts a flux ratio of [CII] to CI of ~ 600, and the in- 
tensity of [CII] and CI would be much higher than the CO inten- 
sities. We note, however, that the fine-structure lines are already 
optically thick at these densities, while the CO lines are optically 
thin. Thus, additional clouds may be present that enhance the CO 
intensities, while the optically thick emission in [CII] and CI is 
unaffected. Nevertheless, as observations show a corresponding 
line ratio of 17.7, these lines should originate from a different 
gas component. Our models require the presence of additional 
gas clouds with density of ~ 10"* cm""*. Go ~ 10''^^ and large 
columns Nh > 10^' cm"^, at which enough of the soft-UV flux 
is shielded to yield a low ratio between [CII] and CI. This line 
ratio is indeed strikingly low, as PDR models generally predict 
larger numbers for this ratio, and also in this case the match is 
not perfect, but only within the 20% error of the CI measure- 
ment. Although this may be an issue due to uncertain abundance 
ratios, a more precise determination of the flux in this line is 
highly desirable. 

So far, we have postulated two gas components to explain 
the CO line SED and the fluxes from [CII] and CI. It is however 
important to cross-check whether these components may affect 
the line ratios that the other component should reproduce. For 
the CO line SED, indeed the intensity in a low-density cloud 
is smaller by one order of magnitude. The high-density compo- 
nent, however, gives rise to an intensity in [CII] that is smaller 
than the corresponding intensity from the other component by 
just a factor of 5. To avoid that this perturbes the [CII] to CI flux 
ratio, we need to require that the low-density gas is more abun- 
dant than the gas at high densities. In this case, it seems however 
likely that the low-density gas would perturb the CO line SED 
significantly. 

The most probable explanation in terms of a two-component 
picture is thus that indeed the low-density component is more 
abundant and explains the observed fine-structure lines. The im- 
pinging soft-UV flux must be moderate in order to reconcile the 
low ratio between these lines. To still explain the observed CO 
line SED, we need to require that the high-density component is 
at higher temperatures, due to an increased value of Go ~ 10^. 
As one generally expects that high-density gas is exposed to a 
weaker radiation field due to shielding eff'ects, this component 
should be spatially separated in regions of strong active star for- 
mation. Alternative scenarios are however feasible as well. For 
instance, the observed [CII]/CI ratio may also be produced in 
the presence of a weak X-ray background rather than a soft-UV 
field, or if a the cosmic -ray background is enh anced by a factor 
of 10 -100 compared to the Milky Way (see also lMeiierink et aU 
I2006 h. 

Of course, this model is still an oversimplification, as one 
may indeed expect a variety of different densities in both star- 
forming and more quiescent parts of the galaxy. In addition. 



there are theoretical uncertainties concerning the metallicity and 
the abundance ratios that may affect our results. Nevertheless, 
this model reproduces the main features in the observed fluxes 
and illustrates the difficulties in simultaneously modeling obser- 
vations in different lines. This challenge will become more se- 
vere, as the high sensitivity of ALMA will allow us to study 
even more lines. On the other hand, it is certainly desirable if 
some of the simplest models can be discarded on such grounds. 
We note that such models are also predictive. For the scenario 
described here, we would for instance expect that also the flux 
in the [OI] 63 /im line and the [OI] 145 fim line originates pre- 
dominantly from the low-density gas component. In this case, 
the ratio of these fluxes to the flux in [CII] would be 0.06 - 0.1 
and 8 x 10"^, respectively. If the oxygen abundance is enhanced 
by a factor of 4 compared to the galactic abundance, the flux in 
the 63 fim line would be larger by a factor of 2. 

With future ALMA measurements, we expect that the CO 
line SED can be probed in more detail and with higher accu- 
racy at higher-7 levels, to check whether the turnover has already 
been reached. Apart from the observed gas component that in- 
dicates a turnover near the J - 6 transition, a higher density 
component may be present that can more easily excite flux at 
higher-7 levels. The increased sensitivity and angular resolution 
of ALMA may help to detect spatial variations in the different 
fluxes, which may help to discriminate regions of intense star 
formation from more quiescent zones. In the center of the galaxy, 
ALMA can check for the presence of an X-ray dominated re- 
gion. Significant progress may however be possible in the mean 
time, for instance by measuring additional high-/ CO lines or by 
a more accurate determination of the [CII] to CI flux ratio. 

6. The expected number of sources 

As shown in § |3] large-scale PDR fluxes, for instance in the 
[CII] 158 fim line or the [OI] 146 fim line, are sufficiently large 
for detection in a few hours. Fluxes from the XDR originate from 
a smaller region, but may also give rise to a significant flux com- 
ponent. With an integration time of a few hours, active galaxies 
with a comparable brightness as NGC 1068 should be detectable 
in an ALMA field of view. One might therefore wonder whether 
a search for active galaxies, based on the before-mentioned fine- 
structure lines, is feasible. Beyond that, it is of strong interest to 
know whether there are reasonable chances to find active galax- 
ies with other surveys such as JWST. In such cases, ALMA 
could perform relevant follow-up studies that probe the central 
X-ray dominated regions and the gas dynamics within them. If 
a sufficient number of sources is obtained, ALMA can probe the 
feeding of black holes between redshift 6 and 10 and thus con- 
strain models concerning the growth of supermassive black holes 
based on the diagnostics provided here. We therefore conclude 
this paper with an estimate on the high-redshift black hole pop- 
ulation. 

To assess this possibility in more detail, we start with some 
general considerations regarding the black hole population at z > 
6. Then we discuss estimates on the number of sources and the 
possibility to detect them with an ALMA deep field. 

6.1. Black hole growth at high redshift 

As shown recently by IShankar et alj (|2009|) . a supermassive 
black hole with 10^ Mq today had, on average, 10^ Mq at z = 5. 
Therefore, one needs to explain how supermassive black holes 
have accreted this mass at early times. Indeed, some of them 
need to accrete even faster, as 10^ Mq black holes have already 
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Fig. 16. The average accretion history of a black hole with 
10^ Mq at z = 5, depending on a coefficient p which is a func- 
tion of the conversion of rest-mass energy to luminous energy, 
the Eddington ratio and the duty cycle. 



been detected at z > 6 (e.g. iFan et all l2006l) . We will however 
focus on the more conservative case of 10^ Mq black holes. 

We follow Shapiro (2005^ and describe black hole growth by 
accretion using the formula, 



dMn 



Mr 



dt 



where 



^growth 



(2) 



(fM/O-l) 1 1 

I - Cm /duty 

s 0.0394/7 Gyr, 



Tgrowth = 0.0394 



Gyr 



(3) 



and where - L/Mqc^ is the efficiency of conversion of rest- 
mass energy to luminous energy, Ae the Eddington ratio and /duty 
the duty cycle, i. e. the fraction of quasars active at a time. The 
parameter p summarizes the complicated depende nce on the lat- 
ter p arameters. With the parameters inferred by 'Shankar et al.' 
(1200 9) for the observed high-redshift black hole population up 
to z = 6, we obtain p ~ 25. 

The dependence of the black hole mass evolution in this pa- 
rameter is given in Fig. [16] For p > 7, the black hole mass hardly 
evolves with redshift, implying that black holes would need to 
form with extremely high masses. There are, however, no theo- 
retical models available that would explain such massive seeds, 
which would have comparable masses to the first galaxies at high 
redshift. 

It is more likely that the black hole population evolved with 
redshift, implying that a larger fraction of them was active and 
had higher Eddington ratios. For instance, a value of p - 2.5 
requires intermediate mass black holes with 10"* - 10^ Mq 
to be the progenitors of the first supermassi ve black holes. 
Such sc enarios have been sugge ste d b y Eisenstein & Loeb 
(Il995h. iKoushiappas etal] (|2004. Ifiegelman et ah (i2006) . 



Spaans & Silkl OOOeh . iLodato & NataraianI (l2006l) and 
Dijkstra et ST (120081). For even smaller values of p, su- 
permassive black holes could even originate from stellar 
progenitors. 

In this respect, it is speculated that so-called dark stars, 
which would be stars p owered by dark mat t er annihilation rathe r 
than nuclear fusion (Spolv ar et all l2008t iFreese et all l2008h . 
could obtain even larger masses than conventional Pop. Ill stars. 
However, the high masses suggested in these works are not 
confirmed by other authors (see also work by lloccol l2008t 



Fig. 17. Estimates for the number of sources in a redshift interval 
of Az = 0.5 with in a solid angle of (T) ^. W e show extrapo lations 
of the results by IShankar et all (120091) and iTreister et all (|2009|) 
to high redshift, for luminosities larger than 10** erg s~'. In ad- 
dition, we give an estimate based on the number of high-redshift 
black holes required in order to produce the present-day 10^ M© 
black holes. As discussed in the text, the actual number of black 
holes should lie in between these estimates. 



Ilocco et"aL l l2008h and m ay be in conflict with the o bserved 
reionization optical depth dSchleicher et alll2008ll2009l) . 

In all of these cases, it is evident that the black hole would 
need to grow considerably by accretion, and that a change in 
the parameters describing the black hole population is required. 
In particular, a larger fraction of active quasars and a higher 
Eddington efficiency are required to obtain th e black hole masses 
that w e observe today. As pointed out by iKawakatu & Wadal 
(119991) . even super-Eddington accretion may be required to ob- 
tain the observed black hole masses. 

6.2. Estimates on the number of sources 

To estimate the expected number of high-redshift black holes, 
we will show extrapolations of theoretical models describing the 
observed high-redshift black hole population, as well as esti- 
mates based on the number of 10^ Mq black holes in the local 
universe. 

To translate a given black hole number density into the num- 
ber of sou rces in a certain field of view , we adopt a similar for- 
malism as lChoudhurv & Ferraral (12007 ) and count the number of 
black holes N{z, Mbh) with masses larger than Mbh in a redshift 
interval [z, z + Az] per solid angle. This is given as 



r'+^' dv r*'"'" 

A^(z,Mbh) = dz'—— flflogM^j, 

J, dz dLl Jm^„ 

X fa- 



dn 



'<i log Mbh' 



(4) 



where dV dz' ' dil ' denotes the comoving volume ele ment per 
unit r edshift per unit solid angle, which is given as dPeeblesL 
[T99I 



dV 
dz'dD. 



Dic- 



dt 

dz 



(5) 



In this expression. Da is the angular diameter distance, c is the 
speed of light and dtjdz - l/(//(z)(l +z)), where H(z) is the 
expansion rate as a function of redshift. The term dn/d log Mi, 
describes the number density of dark matter halos per unit mass, 
and the maximum black hole mass is taken as Mmax - lO** Mq. 
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As explained above, we focus here on the possibility to 
search for fine-structure line emission. With a bandwith of 
8 GHz, we expect that ALMA could scan a redshift interval of 
6z - 0.2 at z ~ 8 in bands 6 and 7 within an integration time 
of a few hours. We assume here that this process is iterated until 
a total redshift interval of Az - 0.5 is covered, which we adopt 
as a fiducial redshift range in our calculation. Alternatively, one 
could focus on the detection of continuum flux, which would al- 
low to scan an even larger redshift interval within a given time, 
but at the same time it would be more difficult to infer the red- 
shift of source. 

For our first estimate th at extrapolate s theoretical models 
of lShankaretaLl |2009.) and iTreister et alj (.2009) , we use their 
results for quasars with luminosities larger than lO^'* erg s 
which can be associated with a black hole mass of 10^ Mq 
for an Eddingt on ratio of 10%, as they find in their calcula- 
tions. Note that lTreister et alJ (l2009l) only provide numbers for 
Compton-thick quasars, but give a Compton-thick quasar frac- 
tion of ~ 10%, which we use to infer the total population. The 
results based on this approach are given in Fig. [17] and are only 
approximate. 

The other approach is based on the number of 10^ Mq black 
holes in the local universe. As noted ab ove, such black hole s 
had an average mass of 10^ Mo at z = 5 (IShankar etalil2009h . 
Therefore, the comoving number density of such black holes can 
be estimated using the number de nsity of today' s lO^M© black 
holes, which is ~ 10"^^ Mpc"^ (IShankar et all l2009l) . It may 
be even high er, as Cornpton-t hick black holes at z = are hard 
to detect (Tre ister et alll2009l) . The expected number of sources 
based on this estimate is shown in Fig.fTT] 

One may note that the number of sources obtained in this 
way is larger th an the value from the extr apolation of theoretica l 
models, i. e. bv IShankar et"an (120091) and ITreister et a n (120091) . 
Such a behavior can be expected from our discussion in § 16.1! 
As supermassive black holes should have hi gher duty cycles at 
z > 6 , the models by IShankar et aH (l2009l) and ITreister et aP 
(|2009|) will naturally underestimate this population. At the same 
time, the estimate based on the present-day 10^ Mq population is 
also an upper limit, since not all of these black holes will be ac- 
tive at the same time. We therefore expect that these approaches 
yield an upper and a lower limit and that the actual black hole 
population will lie in between. 

For the expected abundance of 10^ Mq black holes at z = 5, 
we therefore adopt a fiducial number of ~ 0.08 in a solid an- 
gle of (1')^, and a number of ~ 0.03 at z = 8. We estimate 
the uncertainty to be half an order of magnitude. Based on the 
large fluxes found in § 15. H and the results of Spaans & Meiierink 
(12008). we argue that even smaller syste ms with 10^ Mp black 
holes should be detectable. Calculations bv lShankar et al.l (l2009l) 
show that the cumulative space density of high-redshift quasars 
increases by almost an order of magnitude if the minimum lumi- 
nosity is lowered by an order of magnitude. Additional sources 
may be present as well, such as starburst galaxies or ULIRGS. 
One may therefore expect one active galaxy in a field of view of 
one arcmin^ at z ~ 6. 

At higher redshift, this number may however decrease signif- 
icantly. Therefore, searching for new sources in an ALMA deep 
field seems difficult for z > 6, and it is important that surveys 
as performed by JWS10 provide a catalogue of high-redshift 
sources. Based on the estimates discussed here, we conclude that 
this may be sufficient to study black hole accretion between red- 
shifts 6 and 10, and perhaps even beyond. 

* http://www.jwst.nasa.gov/ 



7. Conclusions 

In this paper, we have explored how ALMA observations may 
help to extend our current knowledge concerning the observa- 
tions of high-redshift quasars. An obvious result is that current 
observations of the starburst component can be done with higher 
resolution and higher sensitivity, which may therefore provide 
detailed local gas kinematics. Due to the higher sensitivity, it 
seems likely that also a number of the weaker fine-structure lines 
may be detected, like [Nil] 122 for z > 2.5, or [SIII] 34 jum 
and [Sill] 35 for z > 11, providing a valuable independent 
probe of the gas chemistry. 

Due to the combination of high resolution and high sensitiv- 
ity, ALMA may however go beyond that and try to explore the 
centers of quasar host galaxies at high redshift. As the models 
presented here indicate, one may expect that the chemistry in 
these regions is dominated by X-ray emission from the central 
engine. In contrast to soft-UV photons. X-rays have high heating 
efficiencies and low efficiencies for dissociation. They can there- 
fore give rise to an entirely different molecular cloud chemistry, 
and can excite high-/ CO lines well above the (10-9) transition. 

Additional probes for the central XDRs are available as well, 
in particular fine-structure lines like [OI] 63 jjm, [OI] 146 fim 
and [CII] 158 //m. In particular for high X-ray luminosities, the 
fine-structure lines of [OI] 63 /im and [CII] 158 jum may provide 
a significant amount of flux even for gas densities of 10"* cm"-'. 
In general, the fine-structure lines become optically thick more 
easily than the rotational CO lines, in particular at densities of 
10^ cm"^. Therefore, they are almost insensitive to the column 
density, but quite sensitive to the X-ray flux. The CO lines, on 
the other hand, have a more non-trivial dependence on the X- 
ray flux and the column density, as in particular for high fluxes 
some dissociation effects will be present. If ideally both the CO 
lines and the fine-structure lines are observed, it may thus be 
possible to derive constraints both on column densities and the 
local X-ray flux. The line ratios may also provide valuable infor- 
mation on cloud gas d ensity and temperature, as pointed out by 
iMeiierinketaP (120071) . 

We have applied our models to available data for observed 
XDRs. For galaxies such as NGC 1068, we expect that the 
CO line intensity rises continuously up to the (17-16) tran- 
sition. Such high ex citation is not possible in PDR models 
("Sp aans & MeiierinH 12008). A detection of such lines therefore 
provides a unique diagnostic for an XDR. 

We consider another example, the lensed galaxy 
APM 08279. There, the situation is not fully clear and is 
currently consistent with the presence of an XDR, which should 
also be expected due to the high X-ray luminosity, but also with 
the presence of a PDR with a radiation field of Go = 10"*^. To 
some degree, this is due to uncertainties in the measured line 
fluxes. We expect that improved measurements with ALMA 
may help to distinguish these scenarios further and probe the 
local gas motions in more detail. 

We finally discussed the chances to find additional sources 
with ALMA in a survey that scans a few arcmin^. General con- 
siderations regarding the average accretion history indicate that 
quasar duty cycles should be larger at high redshift, so that black 
holes can grow to the required masses at z = 6. We therefore 
expect more sources compared to what one expects from a naive 
extrapolation of the observed quasar population. Our estimates 
indicate that about one source per arcmin^ may be present near 
z ~ 6. At larger redshift, the number of sources may decrease 
significantly and one may need to rely on follow-up observa- 
tions of sources detected by JWST. Due to the potentially high 
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spatial resolution, the transport of gas to the center of the host 
galaxy may be probed for the first time in such galaxies even at 
z~ 10. 

We finally note that the diagnostics considered in this paper, 
in particular the high-7 CO fines, may not only be important 
for high-redshift observations with ALMA, but also for future 
space-borne observations in the local universe with SPICA or 
FIRI. 
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